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;_( ' Abstract 

We study Bianchi type / cosmological model in the presence of 
magnetized anisotropic dark energy. The energy-momentum tensor 
^ ■ consists of anisotropic fluid with anisotropic EoS p = ujp and a uniform 

O , magnetic field of energy density ps- We obtain exact solutions to the 

field equations using the condition that expansion is proportional to 
T^ ' the shear scalar. The physical behavior of the model is discussed with 

ly-^ . and without magnetic field. We conclude that universe model as well 

O I as anisotropic fluid do not approach isotropy through the evolution of 

the universe. 
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1 Introduction 

Recent cosmological observations contradict the matter dominated universe 
with decelerating expansion indicating that our universe experiences acceler- 
ated expansion. The accelerating expansion of the universe is driven by mys- 
terious energy with negative pressure known as Dark Energy (DE). The evi- 
dence of the existence of DE comes from the Supernova observations [H E] and 
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other observations such as cosmic microwave background (CMB) anisotropics 
measured with WMAP satelhte [5] and large scale structure |1]. These ob- 
servations suggest that nearly two-third of our universe consists of DE and 
the remaining consists of relativistic dark matter and baryons ^J. 

In spite of all the observational evidences, the nature of DE is still a 
challenging problem in theoretical physics. A variety of possible solutions 
such as cosmological constant [B], quintessence |7], phantom field [S], tachyon 
field [9], quintom [10], and the interacting DE models like Chaplygin gas [TT] . 
holographic models [12] and braneworld models [13] etc. have been proposed 
to interpret accelerating universe. However, none of these models can be 
regarded as being entirely convincing so far. 

Recently, many authors have studied the Bianchi type / model in the 
presence of anisotropic DE. Rodrigues [H] constructed a Bianchi type / 
ACDM cosmological model whose DE component preserves non-dynamical 
character but yields anisotropic vacuum pressure. Koivisto and Mota [T5| [T6] 
proposed a different approach to resolve the CMB anisotropy problem; the 
earlier isotropy of the universe could be distorted by the direction depen- 
dent acceleration of the later universe. Koivisto and Mota [TB] investigated 
the Bianchi I cosmological model containing interacting DE fluid with non- 
dynamical anisotropic EoS and perfect fluid component. They suggested 
that if the EoS is anisotropic, the expansion rate of the universe becomes 
direction dependent at late times and cosmological models with anisotropic 
EoS can explain some of the observed anomalies in CMB. 

Mota et al. [T^ explored the possibility of using the cosmological obser- 
vation to probe and constrain an imperfect DE fluid. They concluded that 
a perfect fluid representation of DE might ultimately turn out to be a phe- 
nomenologically sufficient description of all the observational consequences 
of DE. However, one cannot exclude the possibility of imperfectness in DE. 
Akarsu and Kilinc [TH [19] suggested that anisotropic fluid must not neces- 
sarily promote anisotropy in the expansion whereas such fluid may also act 
to support isotropic behavior of the universe. It has been shown [18j that 
anisotropic Bianchi I model in the presence of perfect fluid and minimally 
interacting DE shows isotropic behavior for the earlier times of the universe. 

Primordial magnetic fields can have a significant impact on the CMB 
anisotropy depending on the direction of field lines [201 121]. Many people 
have investigated the infiuence of magnetic field on the dynamics of uni- 
verse by analyzing anisotropic Bianchi models. Milaneschi and Fabbri [22] 
studied the anisotropy and polarization properties of CMB radiation in ho- 



mogeneous Bianchi / cosmological model. Jacobs [23] explored the effects of 
a uniform, primordial magnetic field on Bianchi type / cosmological model. 
He concluded that the primordial magnetic field produced large expansion 
anisotropies during the radiation-dominated phase but it had neghgible ef- 
fect during the dust-dominated phase. King and Coles [21] discussed the 
dynamics of magnetized axisymmetric Bianchi / universe with vacuum en- 
ergy. He examined the behavior of scale factors perpendicular and parallel 
to the field lines. Roy et al. [23] investigated Bianchi type / cosmological 
models containing perfect fluid and magnetic field directed along x axis. Ex- 
act solutions are obtained using the condition that expansion is proportional 
to shear scalar. 

In this paper, we would like to investigate the effects of magnetic field on 
the dynamics of anisotropic Bianchi / model in the presence of anisotropic 
DE. The paper has the following format. In section 2, we present anisotropic 
Bianchi type / model and formulate the dynamical field equations which 
describe the evolution of the universe. In section 3, we obtain exact solution 
to the field equations and discuss the physical properties of the solution. 
Section 4 contains a brief discussion related to two special cases for /3 = 
and m = 1. Finally, in section 5, we summarize the results. 

2 Bianchi / Model and the Field Equations 

The line element for the spatially homogeneous, anisotropic and LRS Bianchi 
type I spacetime is given by 

ds^ = de - A^{t)dx^ - B\t){dy^ + dz^), (1) 

where the scale factors A and B are functions of cosmic time t only. For 
A(t) = B{t) = a(t), this reduces to the fiat FRW spacetime. This spacetime 
has one transverse direction x and two equivalent longitudinal directions y 
and z. We assume that the universe is filled with anisotropic fluid, and that 
there is no electric field while the magnetic field is oriented along z axis. The 
scale factor A{t) is in the transverse direction, perpendicular to magnetic 
field while B{t) is along the direction of field lines. King and Coles p^ and 
Jacobs [23] used the magnetized perfect fluid energy-momentum tensor to 
discuss the effects of magnetic field on the evolution of the universe. 

Here we take a more general energy-momentum tensor for the magnetized 



anisotropic DE fluid in the following form 

T^ = diag[p + ps, -Px + Pb, -Py - Pb, -Pz - Pb], (2) 

where p is the energy density of the fluid; px, Py and pz are pressures on 
x, y and z axes respectively and ps stands for energy density of magnetic 
field. The anisotropic fluid is characterized by the EoS p = cop, where u is 
not necessarily constant [22]. From Eq.Q, we have 
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T^ = diag[p + ps, -{(^ + S)p + ps, -(w + -f)p - ps, -(w + -f)p - ps] 

where cux = <^ + S, Wy = w + 7, and w^ = w + 7 are the directional EoS 
parameters on x, y and z axes respectively. 6 and 7 are the deviations from 
a; on X and y, z axes respectively. If the deviation parameters are zero, 
then Eq.([2]) represents the energy-momentum tensor for the isotropic fluid 
and magnetic field [5T]. For zero magnetic field, Eq.(I5|) is reduced to the 
energy- momentum tensor of anisotropic fluid [18] . 
The Einstein field equations are given by 

where R^j^y is the Ricci tensor, R is the Ricci scalar and T^^^ is the energy- 
momentum tensor for magnetized anisotropic fluid. For the Bianchi type / 
spacetime, the field equation take the form 

AB 52 

2^ + ^ = P + PB. (5) 

B B^ 

2- + — = -{uj + 5)p + pB, (6) 

A B AB . . .^. 

I+5 + I5 = -i^ + l)P-PB. (7) 

where dot denotes derivative with respect to time. The energy conserva- 
tion equation, T^^^ = 0, leads to two equations for the anisotropic fluid and 
magnetic field [21] 



p+(l+a;)p(4 + 2|)+p(5^ + 27-) = 0, (8) 
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+ P(4+2.|) - 


-- 0, 


PB = 


/3 
54 
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The conservation equation for the anisotropic fluid can be decomposed 
into two parts, 

K;, = K;, + rl^;, = 0, (10) 

where rf^ is the last term in Eq. ([8]) which arises due to the anisotropy in the 
fluid and 2^^ represents the deviation free part of the Tj^.^^. Let us take [18] 

<,, = piSj + 27|) = 0. (11) 

Using this assumption in Eq.dH]), we obtain conservation of the perfect fluid 



d + 2^' 

-A B' 



t^^ = p+(l + a;)p(- + 2-) = 0. (12) 



Equation (TTT]) is satisfied either 6{t) and 7(t) are trivially zero or the ratio of 
expansion rate on x axis to the y axis is equal to —2'-f/6. In order to obtain 
more general solution, the deviation parameter on the x axis 6{t) is assumed 
to be [H] 

2BA 5 1 , , 

and hence the deviation parameter on y and z axes is given by 

lAA 5 1 , , 

where 5{t) and 7(t) are dimensionless parameters and n is the real dimen- 
sionless constant that parameterizes the deviation from EoS parameter. The 
anisotropy of the DE is measured using the relation {6{t) — '~^{t))/uj{t) and 
for n = 0, DE is found to be isotropic. 



3 General Parameters and Solution of the Field 
Equations 

Here we define some parameters for the Bianchi / model which are important 
in cosmological observations. The average scale factor and the volume are 
defined as 

a={AB^)l, V = a^ = AB\ (15) 



The anisotropy parameter of the expansion is characterized by the mean and 
directional Hubble parameters and is defined as 
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where 

^ = i(W) = ,„, = i(^ + 2|). 

is the mean Hubble parameter and Hi{i = 1,2,3) represent the directional 
Hubble parameters in the directions of x, y and z axes respectively, and are 
given by 

X A' y B 

The anisotropy of the expansion results in isotropic expansion of the universe 
for A = 0. The physical parameters like scalar expansion B, shear scalar a^ 
are given by 

= < = I + 2|' (17) 

2 3^ B^ ^ ' 

It is mentioned here that any universe model becomes isotropic for the di- 
agonal energy- momentum tensor when t— )'+oo. A— j-O, V— >+oo and 

to° >0 (p>0) [I21ES]. 

In order to solve the field equations, we use a physical condition that 
expansion scalar is proportional to shear scalar. According to Throne [27] . 
observations of velocity red shift relation for extragalactic sources suggest 
that Hubble expansion of the universe is isotropic within about 30% range 
approximately |2Sl El] and red shift studies place the limit ;|^<.30, where a 
is the shear and H is the Hubble constant. Collins [50] discussed the physical 
significance of this condition for perfect fiuid and barotropic EoS in a more 
general case. In many papers [211 |3T]-[33], this condition is proposed to find 
the exact solutions of cosmological models. It is given by 

A = B"", (19) 



where m 7^ 1 is a positive constant. Subtracting Eqs.dS]) and ([7]), it follows 

that 

A AB B B^ _ n A^ AB B^ 2/3 

I^ AB ~ 5~ 52 ~ 3^l2+ AB^B^'^B^' ^ ' 



Applying the condition given in Eq. flT9|) . we have 

2fl+^p'";-^-""j!=- J' , (21) 

3(m-l) B fi3(m-l)' ^ ' 

where / = m'^ + 4m + 4: is a positive constant. Replacing B = f{B), it follows 

that 

#2 2(3m2 - 3 - n/) /^ ^ 4/3 

rf5 3(m-l) 5 53(m-l) ^ ^ 

which has the solution 

where c is a constant of integration. Thus the spacetime reduces to the form 

ds' = A^dB^-B^^it)dx^-B\t)idy^ + dz^), 
dB 

dT^ 



2(3m2-3-ui) fiflT-2 

(Sm^— 3m— ni) 



CT 3(m-l) 6^T 

X {dy^ + dz^). (24) 

where B = T, x = X, y = Y, z = Z . 

3.1 Some Physical Features of the Model 

Now we discuss some physical aspects of the model (l2l|l . The directional and 
the mean Hubble parameters will become 

-2(3m^+3m-6-n0 6/3T^^ i /r, 

H^ = mHy = m{cT ^(^^^ ^— -V'^, 

3m^ — 3m — nl 

H = '!l±l{cT -'''-l^:V-' - /^^" y l\ (25) 

3 ^m^ — Stti — nl 

We see that these quantities are found to be dynamical and take infinitely 
large values at T = for m > 2. The values of Hubble parameters H, H^, Hy 



decrease with the increase in time and approach to zero as T — )■ oo for 
n < 3m +3m-6 ^ The scalc factors are found to be zero at T = and hence 
the model exhibits point type singularity. The volume of the universe model 
is given by 

V = T"+2 (26) 

and anisotropy parameter of the expansion is found to be 

(m + 2)2 ^ > 

The expansion and shear scalar take the form 

e = 3H={m + 2){cT ^(^^ iT^o ^i ^ ' (28) 

(m — 1)2 , -2(3m2+3m-6-ni) 6BT~'^ , , , 

^ = ^ ^ ^ {cT ^^^^ -^r^^ y}- (29) 

6 dr/i^ — 6m — nl 

We note that spatial volume is zero at initial epoch and increases as 
T — > oo. The expansion and shear scalar are infinite at T = and decreases 
with the increase in cosmic time. Thus the universe starts evolving with the 
zero volume at the initial epoch with infinite rate of expansion which slows 
down for the later times of the universe. The component of magnetic field 
reduces the expansion, shear scalar and Hubble parameters. The anisotropy 
parameter of the expansion is found to be constant and becomes zero at 
m = 1 (it will be discussed as a special case). Thus the model does not 
approach to isotropy for the future evolution of the universe. The most 
general form of energy density is found by using Eqs.(|5]) and ( IT6!) as 

P = ^H\l - f ) - ^- (30) 

The energy density of the DE for the model (p^ turns out to be 

-2(3m2+3rn-6-no ^T-^ ('im'^ j^ Qm + Q - ul) , , 

The energy condition p>0 leads to 

-2(3^2 +3^-6-no /3T-^(3m2 + 9m + 6 - n/) 



c(2m, + l)(3?7i2 — 3?7i — nl) 
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Figure 1: Plots of p verses cosmic time t (a) n < +3m-6 ^^^^ ^ ^ 

with varying values of m as follows: solid, m = 2; dashed, m = 3; dotted, 

m = 4. 



This shows that energy density of the anisotropic DE is reduced by the 
magnetic field. This turns out to be infinite at the initial epoch, its value 
decreases with the increase in time and converges to zero as T — )• oo with 
the condition that n < 3m +3m-6 gj^own in Figure 1(a). However, for n > 
3m +3m-6 ^ energy density decreases after big bang but it starts increasing 
and becomes infinite as T — )■ oo shown in Figure 1(b). 

Using Eqs.f l25|) and fl3T]) in Eqs. flT3|) and f lT^ . the deviation parameters 
6(T) and 7(T) become 



5iT) 



1{T) 



2n(m + 2)(cT" 



-2(3m^+3m-6-ni) 
3(m-l) 



3m'^—3m—nl ' 



3(c(2m + l)T" 



3(m-l) 



/3T-4(3m2+9m+6-n0 ^ 



-2{3m'' +3m-6-nl) 

nm{m + 2)(cT 3(m~i) 



3m''^—3m—nl 
3m^—3m—nl ' 



2(3m2+3m-6-ni) 



3(c(2m + l)T ^i(^=^ 



/3T-4(3m2+9m+6-n0 ' 
3m'^—3m—nl 



(32) 



(33) 



The deviation free EoS parameter can be obtained by using the expressions 
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Figure 2: Plots of {6 — 7)/a; verses cosmic time t with n < ^ 

initial epoch (b) for future evolution, with varying values of m as follows: 

sohd, m = 2; dashed, m = 3; dotted, m = 4. 



for directional Hubble parameters and energy density in Eq. (fT2ll 



u{T) 



-1 - 



-2(3m^+3m-6-n0 

' 4/3T-"^(3m^+9m+6-nO _ 2c(l+2m)(3m^+3m-6-n;)T 3(m-i) 

^ 3m,-^— 3m— nZ 3(m— 1) 



(2 + m)(c(l + 2m)T 



3(m-l) 



,9T-4(3m^+9m+6-n;) - 
'im? —Zm—nl 



(34) 



The anisotropy measure of anisotropic fluid is given by 



(5-7) 



2(3m^+3m-6-ni) 



ra(m + 2f{cT 3('"-i) 



6/3T- 



(3m^— 3m— ni) ' 



U 



-2(3m2+3m-6-ni) 



c(l+2m)(3m2+3m-6-2n0 rTn ~"^"";,X,_i7 i 

(m-1) ^ + 



3/3(m-2)(3m2+9m+6-n;)T- 
3m^— 3m— ni 



(35) 



The deviation parameters 5{T) and 7(T) are found to be finite at t = 
and converges to Anim + 2)/(3m^ + 9m + 6 — nl) and —2nm{m + 2)/(3m^ + 
9m + 6 — nl) respectively as T — t- oo and n < ^m -3"^ _ xhe anisotropy 
measure of the DE (5 — '~f)/u is constant at t = shown in Figure 2(a) and 
converges to 2n{2 + m)^/{{2 — m){3'm'^ + dm + 6 — nl)) for the later times 
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of the universe (Figure 2(b)) with the condition that n < ^^ ^ ^"^ and vice 
versa for n > 3m. ~^'" . We note that the anisotropy of DE does not vanish 
throughout the evolution of the universe. Now we check the behavior of u 
for n < 3!2L_z3m Pqj^ ^]-^g earher times of the universe, the deviation free EoS 

parameter of the DE is given by cj = — 1 + ^(m'2+rn-2)' which shows that 
expansion in the universe may be in the quintessence region after big bang. 
Also, w — )■ — 1 + ^^ as T — 7- oo. Hence for later times, u represent EoS of 
cosmological constant for m = 2 and u may result in quintessence region for 
m > 2. Thus the model represents the accelerating expanding universe. 



4 Some Special Cases 

Here we discuss the following two special cases, i.e., model with (3 = and 
model with m = 1. 

4.1 Model with /3 = 

In the absence of magnetic field i.e., /3— j-O, the model i ^^ reduces to the 
following form 

ds' = ZL.. - T'-dx^ - ndy' + dz% (36) 

For this model, the directional, mean Hubble parameters and energy density 
of the DE reduce to the following form 

, —(3m +3m — 6 — riZ) 

H = ^ ^T ^(^=^) , (37) 

o 

-2(3m^+3m-6-rtZ) 

p = c(2m + l)T ^(^^^) . (38) 

We see that dynamical H, H^., Hy and p are infinite for earlier times and 
converge to zero as T — )■ oo provided that n < 3m +3m-6 ^^^^ ^j^g versa. The 
expansion and scalar are found to be 

Q = 3H = ^ ^T ^(^^^) , (39) 

O 

C(fjl — 1)2 -2(3m^+3m-6-n0 

a = -T 3(m-i) , (40) 
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The expansion in the universe is infinite at the initial epoch and decreases 
with the increase in time for n < 3m +3m-6 ^ However, if n > 3m +3m-6 
then the universe starts expanding at T = and expands indefinitely as 
T — )■ oo. The anisotropy parameter of expansion and anisotropy measure of 
fiuid become 

A = 2^ '-, 41 

(m + 2)2' ^ ' 

(5-7) _ c(m-l)(m + 2)3 ^^^^ 



ijj 



(2m + l)(3m2 + 3m - 6 - 2ln) ' 



We see that the anisotropy of the expansion and that of DE are constant so 
that both model and DE fiuid remain anisotropic. The deviation free EoS 
parameter is found to be 1 — (2Zn)/(3(— 2 + m + m^)) which is a constant. 
bj may begin in the quintessence or phantom region depending on the values 
of constants. 

4.2 Model with m = 1 

For this value of m^ Eq.( TT9l) implies that 

Ait) = B{t) = a{t) (43) 

and the spacetime (II]) becomes 

ds"^ = dt^ - a^{t){dx^ + dy"^ + dz^) (44) 

which is the spatially fiat FRW metric representing homogeneous and isotropic 
universe. Using Eq.( H3|) in (!20|) . we have 

which yields 

a{t) = V2U l^t + cY'\ (46) 

The Hubble parameter and expansion scalar for this model become 

H =-= ^ (471 

a {2l3t + cy/Mj3)' ^ ' 

Q = 3H= ^^-^=r (48) 
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while the energy density is 

3(3(2 -n) , , 

O = — : U9) 

^ 8(3t^ + 12nc'^ + 8cty/6np' ^ ' 

We find that Hubble parameter and expansion scalar are constant at 
initial epoch i.e, at t = and decrease with the increase in time. Hence the 
universe is expanding for the earlier times of the universe. The energy density 
is also constant for the earlier times of the universe and decreases with the 
increase in time. The deviation free EoS parameter is 1/3 which represents 
the radiation dominated phase of the early universe [231 El] • Equation f H^ 
implies that for the radiation dominated universe a(t)cxt^/^ [551 EB], the model 
represents expanding universe. The anisotropy parameter of the expansion is 
zero since Eq.( H^ is an isotropic model. The deviation parameters 5(t), 7(t) 
are constants and hence the anisotropy of the DE —^ = j^z^ is also constant. 
The anisotropy of the DE vanishes if we choose the dimensionless constant 
n to be zero. 



5 Summary and Conclusion 

In this paper we have constructed Bianchi I cosmological model with magne- 
tized anisotropic DE fiuid having anisotropic EoS. The deviation parameters 
6 and 7 are obtained by assuming that conservation equation of DE consists 
of two separate conserved parts. The exact solution of the field equations is 
obtained using the condition that expansion scalar is proportional to a. 
We have discussed some physical aspects of the model both in the presence 
and absence of magnetic filed. 

The component of magnetic field reduces the energy density, expansion, 
shear scalar and Hubble parameters. The expansion in the universe is found 
to be infinite at the initial epoch which decreases with the increase in time. In 
the absence of magnetic field, a similar behavior of expansion is observed for 
n < 3"! +3m-6 ^ however, if n > 3m- +3m-6 i^j^gj^ ^j^g universe starts expanding 
at T = and expands indefinitely as T — )• 00. The anisotropy measure of 
the DE is dynamical and found to be finite for both earlier and later times 
of the universe. The isotropic DE can be recovered by choosing n to be null, 
where n parameterizes the deviation parameters. The universe model does 
not approach to isotropy since anisotropy parameter of expansion A is found 
to be constant with and without magnetic field. For 171 = 1, the Bianchi / 
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model is reduced to spatially flat FRW metric which represents homogeneous 
and isotropic universe for earlier times. For this case, the deviation free EoS 
parameter is found to be 1/3 which represents the radiation dominated phase 
of the early universe [231 El] • 
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